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1.0 INTRODUCTION 

The programming aspects of SFENES are described In the User's Manual. The 
information presented here is provided for the Installation programmer. It is 
sufficient to fully describe the general program logic and required peripheral 
storage. All element generated data is stored externally to reduce required 
memory allocation. A separate section is devoted to the description of these 
files thereby permitting the optimization of I/O time through efficient buffer 
description. Individual subroutine descriptions are presented along with the 
complete Fortran source listings. A short description of the major control, 
computation, and I/O phases is included to aid in obtaining an overall 
familiarity with the program's components. Finally, a discussion of the 
suggested overlay structure which allows the program to execute with a 
reasonable amount of memory allocation is presented. 
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2.0 GENERAL PROGRAM LOGIC 

The general organization of the SFENES program is illustrated in Figures 
1 thru 6. The element routines consists of a control program and four 
principal computation phases; namely, 

1. Input Interpretation phase 

2. Spatial integrations and computations 

3. Matrix computations and time assembly 

A. Stress recovery and element output phase 

2.1 Control Program (Figure 1) 

Routine "INTRFC" controls the element computation flow. Initially, the 
"Input Phase" and "Spatial Computational Phase" are executed. These two 
phases are outside of the iterative phases and need only be executed once. 
Routine "INTRFC" then executes the "Matrix Computation and Temporal Assembly 
Phases". Control is referred to the "MAIN" program for all spatial assembly 
and displacement computations, [ Ky ] {u} = {F^} . The "INTRFC" portion is re- 
entered for convergence checking and the stress recovery phase is then exe- 
cuted. The iterative process is repeated until convergence or the maximum 
number of iterations is achieved. The Control program is the interface be- 
tween the element routines and the MAIN program unit. 

2.2 Input Interpretation Phase (Figure 2) 

The purpose of this phase is to read the formatted element input data and 
control data from the input file "IOIN". The data, in some cases (element 
loads and local grid coordinates), is used in additional computations and is 
finally output to binary files "I0EL", "IOLD", and "IOUS" for use by the 
remaining phases. This phase is only entered prior to the first iteration of 
the problem. Routine "INPUT' reads the formatted data, provides an echo print 
for the user and outputs all data to the output files. "GRIDG" computes all 
of the remaining grid point coordinates from local grid points 13-16 and the 
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corner thicknesses. Routine "LOADIN' 1 reads the input element load data and 
computes the resulting grid point forces at the user input time increments. 

2.3 Spatial Integrations and Computations (Figure 3 6 3A) 

The "Spatial Computation Phase" is executed once prior to the first iter- 
ation* The function of this phase is to perform all spatial computations and 
any other "one" time computations which required spatial integrations* The 
primary tasks of this phase are outlined briefly below. 

a) The volume, area, and surface definite integrals are computed for the 
polynomial terms xPy®z^ and p^x^y^z^; 0<n<9, 0<m<9, 0<1<7, 0<k<3, 0<i<7, 
0<1<7. These integrals are computed by routines VOLCOF, COLXYZ, CEAXYZ, 
CNLXYZ , CURCOF, CURPSI, CNRPSI, CURXYZ, CNRXYZ, TRLCOF , VOL-, VL-, REA-, 
RA-, SUR-, and SUP-. 

b) Routine LOADD computes the element force vector {F} by combining the 
input element loads as a function of time, (file IOLD) the appropriate 
spatial integrals from (a) above, the Hermitian time functions, and inte- 
grations over the local time ranges. The resulting combination of ele- 
ment body force, pressure and traction forces are output to file "IOFO" 
for intersolution of [Kp]{u}*{F} equation. 

c) The final function of this phase is to compute and store the spatial 
matrix partitions of [M], [K], [X], [A], [B], [ C e ] , [B 2 g]t I p 2 8 ^ 

[N x ], and [ p 28^ These matrices are stored on file "IOCF", 

"ICPR", and "IOTP" for use by the third phase. The matrices computed In 
this task are reused In subsequent iterations since the shape functions 
are constant. 

2.4 Matrix Computation and Time Assembly Phase (Figure 4) 

The function of this phase is to compute the element stiffness and stress 
matrices for each local time range. This involves the computation of the [Eg] 
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matrices as a function of time, (routine EPQP19) and integrating over the 
local time intervals. The result of the time integrations is multiplied by 
the spatial matrices as computed in the spatial computation phase and assembled 
into the element [K e ], [S g ], [A g ] and [B e ] matrices for each local time in- 
terval. These computations may be made for either the linear elastic quad 
plate or the non-linear elastic-plastic quad plate elements (routines LEQP and 
EPQPB respectively). Spatial matrices are input from files "IOPR" and "IOTP" 
and prior iterative values of {w} and {o} are input from file "IOUS". The 
"IOAC/ 1 and "IOAl" are used alternately to read the previous and store the cur- 
rent coefficients of polynomials in (t) used to compute the [Ey] matrices. 
File "IOKE", "IOAE", "IOBE", and "IOLD" are output files for storing element 
[Kg] [A e ], [B e ] and [S g ] matrices respectively. 

Routines "TBASEM" and "TABSEM" perform the temporal assembly of the 
element [Kg] matrices for the "LEQP" anbd "EPQP" elements respectively. Files 
"IOKE", "IOCE", "IOBE", "IOAE" control the input to these routines and the 
final element [K^] matrices are output to file "IOKF". For the "EPQP" element 
the computation [B ]*[C e ]“* is stored on file "IOLD" for element stress 
computations . 

2.5 Stress Recovery and Element Output (Figure 5) 

The "CNVRGV routine checks for convergence of all displacement/velocity 
values and stores the latest value of [u] on file "IOUL" for the next 
iteration. Element stress is recovered in routines "CSIGO" or "CSIGL" for 
element types "EPQP" and "LEQP" respectively. The input to these routines is 
contained in files "IOTF", the current displacement/velocity values {u}, and 
"lOLtf', the stress matrices [S] or [Cg] _ * fB g ]. The displacement/velocity [u] 
and stress {a} vectors are output to file "IOUS" to be used in the computation 
of [Ey] for the next estimate of [Kg]. Printed output of selected grids 
displacement and/or stress values for each iteration is provided for by 
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routine "WRITO". This output is at the major time intervals only. 
Convergence of the problem initiates the extrapolation and output of 
displacement /velocity values at user specified time intervals within the major 
time increments. 

2.6 Displacement Computation (Figure 6) 

Control is returned to the "MAIN" program after the stiffness [Kp] and 
stress ( S ) or [C e ~h*[B e ] matrices have been time assembled in the matrix 
computation phase. The "MAIN" program calls routine "COMPlf' which performs 
the following functions: 

a) Read the partitions of the [ K f ] matrix (see Figure 7 to Figure 9 for 
description of [Kp] matrix) 

b) Applies the user input boundary conditions to the individual partition 
(eliminates rows and columns) and forms the reduced [Kp] matrix. 

c) Eliminates all rows containing displacement related data at times T»0 and 
T=t N (rows 1-84 and rows IN* 168+1 to N* 168+84 of the [Kp] matrix). 

d) Reads the {Fp} vector partitions and time assembles the local partitions 
(see Figure 14 for a description of the assembled {F} vector). 

e) Applies user input boundary conditions as in step b) and c) above to the 
{Fp} vector rows. 

f) Computes {u} = [K F _1 ] {F} and stores data or file "IOTP" for stress 

recovery and output phase • 

The newly computed displacement vector {u} is then read by the stress 
recovery phase and the next iteration is entered if required* 
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3.0 MATRIX DESCRIPTION AND LABELING 

The basic matrices generated by the element computations consists of the 
[Kp], [S], and {Fp} matrices. The matrices are stored on files "IOKF", 
MpOLD'^ and "IOFO" respectively. The basic matrices are illustrated in 
Figures 7 to 14 and are described as follows. 

3.1 Stiffness [Kp] Matrices 

Figure 7 depicts the time assembled fKp] matrix for the EPQP element. 
Each numbered box represents a 168 by 168 partition and each partition is 
labeled from T Q to T N (global time, were Tq-0.0). The numbering of the 
partitions signifies the order in which they are stored on the sequential file 
"IOKF”. Similarly, Figure 8 depicts the time assembled [Kp] matrix for the 
LEQP element. It should be noted that rows 1-84 of partition rows Tq and T^ 
are to be reduced from the matrix before using to solve for {u}. Finally, 
Figure 9 illustrates the labeling of each Kp partition. 
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3.1.1 [K f ] matrix labeling 

The [ Kp ] matrix partitions are computed and time assembled (Kp «K g + 
[ A e ] [C e -1 ] [B e ]). The non-zero 168 by 168 K F partitions are stored on File 
IOKE. Figure 7 shows the Kp matrix-partition labeling and order of storing on 
file IOKF for the EPQP element. Similarly, Figure 8 represent the Kp matrix 
for the LEQP element. Figure 9 is a sketch of the labeling of each partition 
for both the EPQP and LEQP elements. 


FIGURE 7: EPQP [Kp] MATRIX PARTITIONS 



where : 


Note : 


T^l) = input global times 
T q = 0.0 

1) Matrix partition col. no. 1 contains time Tj partitions 

2) Partitions are output to file IOKF in the order they are 
numbered in diagram 

3) Each partition is a 168 by 168 matrix, see Figure 9 for 
labeling of partitions. 
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FIGURE 8. LEQP [K p ] MATRIX PARTITIONS 
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3.2 Stress [S] Matrices (EPQP Element) 


The equation to compute the six (6) stress components at the 28 element 
grid points for time interval J is: 

{°}j " t c 0 ~ 1 HB 0 1 ]{u 0 } + (1) 

l ([C 0 " 1 ]*([B 0 I ] + [B t ( I - 1 >])*{u ( i _ 1) }) + 

1=2 

( c a 1 J f B T J ] f u ( j)J “ 

l ([c 0 - 1 Hq<i)}> + {o 0 } 

i=i 

The time assembled resulting matrix equation, collecting like terms, 
could be represented as follows: 


\ 
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( N+l ) *168 x 

•168 (N* 168) 


Figure 10 illustrates the [Sj] matrix which is computed in routine 
"TASSEM" and stored on file "IOLD" in 168 by 168 partitions. The number of 
each partition represents its location on the sequential output file "IOLU' as 
well as the terms of equation (1) above as listed below. 
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4, 6, ... (Nx2) 

[ C 0 ^ ] [B j 

5, 7, ... ( Nx2)-1 

[C 0 _1 ]*([B ( 
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Figure 11 illustrates the [ S 2 ] matrix. Matrix [ S 2 ] is simply the time 
assembled -[C 1 ] -1 matrix where with the first column removed; i.e., 



File " IOCE" contains a 168 x 168 partition, [C 0 _1 ]. The remaining 
variables namely, {w}, { q } , and { u o / 0 o ) are stored on files "IOTP", "IOQT", 
and "IOUS" respectively. 

Figure 12 contains a sketch of the partition labeling for both the [Sj] 
and [ S 2 ] matrix partitions. 
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FIGURE 11. STRESS [S 2 J MATRIX 
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NOTE: 

1) Each partition is identical and the 168 x 168 matrix [C^ *] is stored on 
file " IOCE" by routine "EPQPA" . 


21 






3.3 Stress [S] Matrices (LEQP Element) 

The 6 x 84 x 28 stress matrices for the LEQP element are generated in 
routine LEQP and stored on file "IOLD". The equation used to obtain [S] is 
[SJ*[E 6 ] [N x J ](grid)’ where I^ 1 ] is evaluated at each of the 28 grid points of 
the quadrilateral plate. Figure 13 illustrates the labeling of [Sj. Since 
[S] for the LEQP element is time independent, the same matrix is used for each 
time increment in the problem solution and the velocity terms of the {u} 
vector are not used in this solution. 


FIGURE 13. LEQP STRESS MATRIX LABELING 
FOR ALL TIME INCREMENTS 
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3.4 Element Force Vector (Fy) 

Figure 14 illustrates the force vector labeling and the time assembly of 
the element load vectors. The equations to adjust the force vectors for known 
initial conditions are as follows: 

< F F(o)> “ O - ^oo 1 ) < u o> - <0 <°o> (1) 

{F f ( 1)} - {F q 2 + F^. 1 } - ( [ K To 1 ] + [A 0 2 ] [C 0 _1 ] [B 0 1 ] ) {u 0 } - (2) 

([A q 2 ] + [At 1 ]) {q 0 } + [a 0 2 ][c 0 " 1 ] (q j } 

If N=1 {F f(1) } - {F f(1) } - [At 1 ] [C a _1 ] [B 0 1 ] {u 0 } 

{F F (j)} = {F 0 (J+1) +F t J } - ([A 0 (J+1) ] + [At J ]) 1C-') [B d 1 ] {u 0 } - 

J=2,N 

([A q (J+1) ] + [At J ]) P 0 } + 

(([A p (J+1) ] + [A T (J) ]) [C a _1 ] {q L } + 

L= 1 

[A 0 (J+1> ] ICo-'Hqj) 

The first force vector is written to file IOFF^ 1 ' for solution of the 
problem [ K F ] {u} = {F ? }. 


NOTE: 

(1) The above computations are not currently performed in the program. It is 
assumed that {u Q } = {0 0 } - {q} - {0} and therefore file "IOFO” contains 
the {F oit L } vectors (unassembled). Routine "COMPU" currently preferring 
the time assembled. 
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4.0 EXTERNAL DATA SET STRUCTURE 


The SFENES program uses 18 data sets during execution. The dataset 
variable unit names are stored in labeled common "!(/' and are initialized to 
the integer unit number in routine "INTRFC". Table 1 contains the dataset 
names, default unit number, routines using the dataset and a brief description 
of the contents of each dataset. Table 2 lists the record descriptions and 
record formats contained in each dataset. All datasets contain binary output 
except "IOIN" and "16". The binary datasets in several instances contain 
variable length records. 
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INPUT OUT 
IPRIN OUT 
PRINI OUT 
WRITD OUT 
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IOKE 10 EPQPB OUT Element [Kg] matrix for each time increment 

(Binary) TASSEM IN 

LEQP OUT 

TBASEM IN 


TABLE 1 (CONTINUED) 
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FILE NAME 
IOLD 

(Load Data) 


(LEQP 

Stress 

Matrix) 


(EPQP 


Stress 

Matrices 


[Cel- 1 * [B e ] 


IOEL 
(Element 
Input Data) 


RECORD DESCRIPTION 

1 

Element Load Code 
2 

Element forces at NGRID 
affected grids for this 
load type and for D.O.F. 
u,v,w 

3 

If ID=2 only element rate 
of force at NGRID affected 
grids for this lead type 
and D.O.F. u,v,w - 


RECORD FORMAT 


( (LCODE(I, J) , 1=1,7) , J=l, 6)' 
42 integers/record 


((FV0(I,J),J«1, NGRID). >1, 3)1 

NGRID = 10, 12, or 28 

.'. 30, 36 or 84 real data/record 


( (FVT (I, J) , J=1,NGRID) , 1=1,3) 

NGRID - 10, 12 or 28 

:30, 36, or 84 real data/record 


NOTE: Records (2) or (2 and 3) are repeated for each input 

load type and for each time point. 


The data is overwritten by: 


1 

Element stress matrix 
6 components, computed 
from 84 displacement 
D.O.F. for all 28 grids 

or by: 


1 

Element stress matrix 
partition [C ]~^([B 0 ]^+ 
[B T ]N-1) 

2 

Element stress matrix 
partition [C ] - l*[Bj]^ 


( ( (STRM(I, J ,K) , 1=1, 6) , J=l, 84) ,j 

K=l,28) J 

14112 real data/record 


( (D0T(I, J) ,1=1,168) , J=l, 168) 
28224 real data/record 


( (DT (I, J) , 1=1,168) , J= 1 ,~168)~~' 
28284 real data/record 


NOTE: Records 1 and 2 are repeated (NTIME-1) times (N=l, 

(NTIME-1) ) 


1 

Element name and code 


ANAM, NELEM 

21 real and 1 integer /record 


2 

Element local coord., 
material prop., and 
load codes 


(X(I) , Z(I)=1,28) , (H(J) , J=l, 4) , 
RHO, AMU, ALPHA, NU, E, YSTR, AJ, 

(ILOAD(T) , 1=1,7) , NONLIN 

94 real and 9 integer data items/ 
record 


TABLE 2 
30 


FILE NAME 

RECORD DESCRIPTION 

RECORD FORMAT 

IOUS 

1 

((SIG(I.J) , 1=1, 28) , J=l, 6) 

(Element 

Element stress components 

168 real data/record 

Input Data) 

(6 each) for 28 grid poin 
points 



2 

((Ul(I,J) ,1*1,28) ,J=1,6) 1 


Element disp-vel. compo- 
nents (u,v,w,u,v,w) for 
28 grid points 

168 real data/record 


NOTE: Records 1 and 2 are repeated NTIME times 

IOIN 

Several 

See User Input Description 

(Element Input 



Formatted) 



I OPR 

1 

((B28(I,J) ,1=1,28) , J=l, 28) ! 

(Spatial 

Computation 

The [ B 28 3 — 1 matrix 

784 real data /record 

Data) 

2 

( ( (DNXPP(I, J,K) ,1=1,7) , J=l, 84) , 

The [ D ] • [ N x ' ] computation 

K=l, 28) 


for all 28 grid points 

16414 real data/record 


3 

((SIGNX(I,J),I=1,28) , J=l,28) 


The {P28^‘l B 28l computa- 

784 real data/record 


tion for 28 grids 

I=polynomial terms 1-28 
J=grids 1-28 


4 

( (P28(I, J) , 1=1,28), J+1,28) 


The evaluation of (P28^ 

784 real data/record 


for 28 grids 

I=polynomial terms 1-28 
J*grids 1-28 


5 

( (CUUM(I, J) , 1=1 , 28) , J=1 , 28) 


The mass computation 

784 real data/record 


[ B 2g -1 ] T '[A/ p 28 T * p 28 dv l ‘ 

[b 28 _1 ] 



6 

(((CIJK(M,N,K,),M=1,28),N=1,28) 


The [K] computation 

K-1,9) 


[ B 2 8 -1 ] T \tAr p 28( L >* 

p I T - p j dV ]*[ B 28 _1 ] 

7056 real data/record 


NOTE: Record 6 is repeated 

for LEQP NLIN=1 
for EPQP NLIN=28 

NLIN TIMES and L=l, NLIN: 
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FILE NAME 


RECORD DESCRIPTION 


RECORD FORMAT 


IOPR 
(Cont . ) 


Record: NLIN + 6 

The [X] computation - 


( i (aBio _1 (j)3 T • 

j=i 


< 


[^p 28 a)* p io T(J) ' p i ds] 

t B 28 _ ^ T " C s P 28 (L) * P 28 T 

P I dS])-ANG<K,J)))-[B 2g 


(J)* 

- 1 ] 


i ( ( (Cl JXE (M, N, K) ,M= 1 , 28 ) , N=1 , 28 ) , I 

|K=1,6) 

4704 real data/record 

K-123456 
I » x y z x y z 
> ANG (K, J) - SIN SIN SIN COS COS COS 


NOTE: Record is repeated NLIN Times and L=l, NLIN; 

For LEQP, NLIN=1 
For EPQP, NLIN=28 


Record: 2*NLIN+6 

The [A e ] spatial computa- 


tions - EPQP only 
f ([B, fi - 1 ] T *[/Pi T -P 2 8dV] + 

T . T 


28 

((B, 


-1 v 
j=l" J 10 (J) J 


'si P 10(J) 


< 


P 28<J) dSI - 'V 1 ' 1 ' 


'su P 28<J) T ' P 28(J) liS]) 
^ ANG (K,J)))* [ b 28 -1 ] 


!( ((CIJAfM.N.K) ,M c l,28) ,N=1,28) ,K=1,3) 
2352 real data/record 

^ K = 1 2 3 

I - x y z 
ANG(K.J) = SIN COS 0.0 




Record: 2*NLIN+7 


( ( (CIJA(M.N.K) ,M=1 , 28) , N=1 , 28,K=I ,3) 


The [B e ] spatial computa- 
tions - EPQP only 


[B 28 1] 


V P 28 (L) ' 


28 


P I dVl * lB 28“ 1] 


K = 1 2 3 

I * x y z 


Note: Record is repeated NLIN times and L=l, NLIN; 

For LEQP NLIN=1 
For EPQP NLIN=28 
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FILE NAME 
IOTP 

Temporary 

Storage) 


IOCE 

([C e -1 ] Matrix) 


IOKE 

([K £ ] Matrix) 


IOAE 

( [ A e ] Matrix) 


IOBE 

( t B e ] Matrix) 


RECORD DESCRIPTION 


Same as File IOPR 
- or - 

Record 1 - EPQP only 
[A q ] n partition of [A e ] 
matrix 

Record 2 - 

[A t ] N partition of [A e ] 
matrix 


RECORD FORMAT 


i((ABIB(I,J). >1.168). >1,168) 
28224 real data/record 


!((ABIN(I,J). 1-1.168). J-l. 168) 
28224 real data/record 


NOTE: Records 1 and 2 are repeated NTIME times 


-or - 


Record 1 

Displacement Vector {u) (U1(K,I),K=1,168) j 

168 real data/record 

NOTE: Record is repeated for 1=2, NTIME 


Record 1 - EPQP Only 
[Ce -1 ] matrix 

Record 1 

[K e ] matrix partition 


;((CE(I, J) . 1=1 . 168) , J=l, 168 
28224 real data/record 

:((AKE(I,T).I-1.336),J=1,336) 
112896 real data/record 


NOTE: Record is repeated (NTIME-1) times 


( (AE(I, J) , 1=1 . 168) . J=1 , 168) 
28224 real data/record 


Record 1 - EPQP Only 
[A 0 ] N partition of [A g ] 
matrix 

Record 2 

[A t ] N partition of [ A^] 

NOTE: Records 1 and 2 are repeated (NTIME-1) times 


( (AE(I.J), 1=169. 336), J-l. 1 68) 
28224 real data/record 


Record 1 - EPQP Only ((BE1(I, J) , 1=1, 168) , J=l, 168) 

[B partition of [B ] 
o £ 

Record 2 ((BE2 (I, J) , 1=1 , 168) , J-l , 168). 

[B ? ] n partition of [B g ] 

NOTE: Records 1 and 2 are repeated (NTIME-1) times 
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FILE NAME 


RECORD DESCRIPTION 


RECORD FORMAT 


IOA1 

IOAO 

(E 7 Coeff 
Arrays) 


IOFO 

(Integ. Force 
Vector) 


IOKF 

([K f ] Matrix) 


IOSE 

([A e ] • [C e ] 
• [B e ] ) 


IOUL 

(u, v, w, 
u, v, w) 


Record 1 - EPQP Only 
Arrays AC, ATI, AT2, AT3, 
AT 4 , All (7, 7, 28, 4) 


(AC, ATI, AT2 , AT3, AT4) ] 

All dimensioned (7, 7, 28, 4) 
27440 real data/record 


NOTE: Record is repeated (NTIMES-1) times 


Record 1 
Vector {f} 


!(( (VL0AD(I, J,K) , 1=1,28) , J=l, 3) , ; 

iK-1,4) ! 

336 real data/record 
1=GRID, J=D.0.F. (u,v,w,) 

K=0, 0, T, T 


NOTE: Record is repeated (NTIME-1) times 


f((ACB(I.J),I=l,168),J=l,168) 
28224 real data record 


Record 1 

Partition (M,N) of the 
[K_] matrix 

NTIME 

NOTE: Record repeated ( E (NTIME - (I-l))-l times 

1=1 


Record 1 - EPQP Only ; ((ACB(I,J),I=1,168)J=1,168) 

Partition (M,N) of 28224 real data/record 

[Aj-tC^l-lBj 

NTIME 

NOTE: Record repeated ( £ (NTIME-I))-1 times 

1=1 

Record 1 (UNEW(I, J) ,1=1,168 

Displacement/velcoity 

vector from the previous 

iteration {u} 


NOTE: Record is repeated for J=2, NTIME 
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SLAVE FINITE ELEMENT FOR 
NONLINEAR ANALYSIS OF ENGINE STRUCTURES 
USER'S MANUAL 


t 


1.0 INTRODUCTION 

The purpose of the SFENES program is to provide the user 
with validation of the basic matrices needed to compute the 
transient displacements of a group of non-linear elements and to 
provide for element stress recovery unpon convergence of the 
problem. In addition, the element matrices routines have been 
provided to compute resultant displacements, checks for conver- 
gence, output displacement and stress and perform element stress 
recovery. The program is limited by the following factors: 

1. The only element currently provided is the quadri- 
lateral plate element with either linear elastic 
analysis or non-linear elastic-plastic analysis. 

2. Only a single structural element may be analyzed 
(multiple time intervals). 

3. All input and output is in terms of the "local" coordi- 
nate system. 

4. A maximum of 25 time points (solution times) may be 
requested. Output displacement times are unlimited. 

The following pages describe the program and its input/ 
output features. 
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2.0 SFENES OVERVIEW 

The Slave Finite Element procedure is described here. 
Reference to Figure 1 will aid in the description of the solution 
algorithm. The element input data is interpreted and stored for 
later reference. All of the required spatial definite integrals 
are computed and stored as is the one time computation of the [C] 
matrix. The input element lead data (if any) is combined and 
integrated over the time intervals. The force vector { F} is then 
stored for later use in computation of { u } . 

The iterative procedure starts with the computation of the 
matrices [K e ] and either [A e ], [B e ] for the EPQP case or [S e ] for 
the LEQP case. For the EPQP case the initial iteration assumes 
that the structural element is elastic for all time intervals but 
for each subsequent iteration the elastic and/or plastic regions 
are computed. The structural matrices are then assembled for 
global time. For the EPQP case final time assembled [Kp] matrix 
is computed as [Kp] * [K e ] = [A c ] [C e ] ^ [B g ] and the matrix 
multiplication [ C e ] ^ [B e ] is stored for later stress 
computation. The linear LEQP case simply involves the time 
assembly of each of the [K e ] matrices as computed over the speci- 
fied local time intervals. The [Kp] matrix and the { F } vectors 
are then reduced based upon the input boundary conditions as well 
as the removal of the tows corresponding to displacements at 
times Tq and T ^ (rows 1-84) and rows ((N*168) + l to (N*168) + 84) 
where N = the number of input time points minus 1 (1<N<24). The 
reduced [Kp] matrix and [F] vector are used to compute the new 
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(START) 


START? 
NEW [Kr] 


(NEW [Kj.]) 


(NEW {u} & [S.]) 


X. INITIALIZE 
2. INTERPRET 
ELEM. INPUT 
AND STORE 


ITERATION 

NO. 


CHECK FOR 
CONVERGENCE 


STORE NEW {u} 
FOR NEXT 
ITERATION 


SPATIAL ASSEMBLY 

1. APPLY BOUNDARY 
CONDITIONS 

2. REMOVE THE T 0 
AND Tji DISPLACE- 
MENT ROW D.O.F. 
[Kf] AND {F} 


COMPUTE: 
ELEMENT FORCE 
(F) VECTOR 
AND STORE 


ELEMENT 

TYPE 


1. COMPUTE , 

{oMcrMBjfiD 

2. STORE 
{0} & {u} 




(NEW [Kj.] AND [S e ]) 
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(u) vector ( [ K F ] {u} = { F } ) . The updated value of { u } is checked 
for convergence and then stored for the next iteration. 

The vector { u } and the stress matrices [C e ]"^ [B e ] or [S e ] 
are used to compute the stress vector {o} * [ C ~ ^ [Be g ] { u } or 
•f a ) = E s e ] {u}. The displacement and stress vectors are then 
stored for use in computation of the [Ey] matrices for the next 
iteration. After each iteration the displacement and stresses 
are printed for the input time intervals only. When the problem 
converges the displacements are also printed at the requested 
output interval. The iterative procedure continues until conver- 
gence or maximum iteration is achieved. 


38 



3.0 QUADRILATERAL PLATE DESCRIPTION 

One finite element is provided in SFENES and that is the 
quadrilateral plate element. Figure 2 depicts the local coordi- 
nate grid labeling and the orientation requirements. The 
analysis procedure permits transient analysis of the plate as 
either a linear elastic element (LEQP) or as a non-linear 
elastic-plastic element (EPQP). 
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4.0 INPUT OVERVIEW 

The input requirements and options are described below. The 

following general rules must be observed. 

1. Each input line consists of 8 fields of 10 characters 
each. In some cases a field is subdivided for integer 
or alpha input. 

2. Integer data and floating point " E" format data must be 
right justified in their respective input fields. 

3. One of the material prop inputs Is required. 

Listed below is a brief summary of the input options and 

notes regarding order and required usage. 

TITLE Title information 

TIME (Required and must preceed and load input) 

The number of time increments +1 and the values of 
beginning/end of the time increments. 

EPQP Contains material properties and local grid point co- 

ordinates for the linear plastic quad plate element. 

LEQP Contains material properties and local grid point co- 

ordinates for the linear elastic quad plate element. 

LOAD Defines element body force, pressure and traction load 

(must be preceeded byTIME input option). 

ISTRES Defines any initial stress component values (at TIME = 

0 . 0 ) 

ISISP Defines any initial displacement and velocity component 
values (at TIME = 0.0) 

PD0F Displacement and STRESS output print selection (Default 

= Print All D.O.F.) 
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PGRD 

ICND 

CONT 

END 


Local grid number output print control (Default = Print 
All Grid Points). 

Boundary conditions, grid identifications (Default 3 
None Bounded Out) 

Execution control parameters - maximum iterations, 
convergence criteria and output format . 

Require last input for the element. 
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5.0 INPUT DESCRIPTION 

Figure 3 contains an illustrative example of typical user 
input. A detailed description of each line and field of input, 
that may be entered for each of the input options summarized in 
Section 4.0, is contained on the following pages. Input record 
description contained on these pages includes: 

1) Input Code, i.e., "TIME" 

2) Record Format, i.e.. Field Usage and Mnemonic 

3) Record Example, as contained on Figure 3 

4) Field Description - verbal explanation of each field in 
the record along with data type, limits and defaults if 

any . 

5) Qualify ing Notes 

All input is formatted and is currently read from logical unit 5 
" IOIN" . 
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0.0005 


0.0010 




)R0 • 16C.0 10.0 0.0 1 

.0 
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• 0 

>RD 1510*0 * 10.0 I 

)A0 , 2 2 1 1 1 1 :v : ■ 

.0 

- 0.0 0 0 0 0 0.0 

k TDCC ' 

0 0 0 0 .-10000.0 0 0 0 

>JKE2fc UZ 

C40.0 C.O 0.0 0 

C2C.C 0.0 0.0 0 

►ISP 03 

• 0 0 • 0 0 *0 

• 0 Q»0 0 #0 


040.0 ' C. 

C20.0- 0. 




OF 

1 1 1 11 1 1 11111 


flj 


11 12 13 14 16 17 16 20 21 


FIGURE 3. INPUT SAMPLE 
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INPUT - TIME (REQUIRED) 


I 

I 


O _J 
3 . <E 

OOlii 

.j a or 


UJ UJ UJ UJ 
T X T X 
hhhh 


^ ^ I I 



I H 1 

I 


1 W 1 

HH | 

Q -J 

1 


o 


I Z 1 

1 


1 O 1 

H 1 

U UJ 

1 
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SO CHARACTERS OF ALPHANUMERIC DATA 


LL O 

0 CD 

LL 0 s 
LU n 


iLi O 
« z 

IL 



I 

l 

f— - 

« 1 



I 

1 



i 

1 

z - 

Z 1 



1 



t 

1 

UJ K* 

Q i 



1 



I 

l 

z 

^ 1 

Q 

-J 

1 



-- -- 

— 

UJ X 

H 1 

-J 

UJ 

1 



1 CM 

I 

-J u 

LL 1 

UJ 

Oi 

1 



1 O 

1 

UJ X 

•-* 1 

f— ( 

<x 

1 

UJ 

J* 

1 

1 

<x 

X 1 

Ll 

-J 

I 

-J 

u 

1 

I 

Ll Z 

U I 



1 

H 

J 

1 

1 

o 

0 1 



1 

>-« 

LL 

1 UJ 

1 

•t 

UJ 1 



1 

H 

z 

1 -J 

1 

J— u 

Q I 



1 


c 

I H 

1 

0 H 

l 



! 


X 

1 *-H 

1 

uj <r 

D 1 

Q 

■ 

1 


LU 

1 H 

1 

i— Q 

.J I 

X 

o 

1 

■H 


— — 

— 


UJ ! 

<X 

z 

1 






| 

u 


1 






Ll 1 



1 



W 

u 

z to 

V' 0 

^ z 


<r 0 

LU i-f 

0 _ J 
UJ 

H H 
Z X 
*-* CL 
~ H 
X 
Z O 
O 

-I > 
-J IL 

o »-* 

Ll H 
Z 
H UJ 
<1 C 
X »-» 
H 

O 
0 H 
= Q 



1 

UJ X UJ 


I 

J <Z X 


1 

HU <E 


l 

M 


I 

H UJ X 


1 

= -1 U 


1 

f— *-* 


1 

UhI 


1 

« H Z 


i 

X 


1 

UJ U. 0 


1 

Z O U 


1 

X 

0 

1 

Z X X 

f- 

1 

<L UJ UJ 

Z 

I 

X 0 Z 

Ul 

1 

X Z X 

H 

1 

X X z 

Z 

1 

c z <r 

o 

1 

X 

u 

1 

UJ UJ CL 


1 

XXX 


1 

H H <1 


0 

I 


U 


w 


c c . 
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U f.l u u u u 


INPUT 


EPDP 


ELASTIC PLASTIC QUADRILATERAL PLATE ELEMENT 


FORMAT s 



EPC-F 
1 _ 

NL 


:e 

!G 

1 r- 

1 AMU 

1 RHO 

“ « 

J ALPHA 

1 A J 

— — # 
1 


! 

YSTR 

NU I 1 1 l l J 


ICORD 

IG 

XENTER 

YENTER 

i ZENTER 

1H 

1 

1 

1 - 

1 


EXAMPLE: 


, * » i "" — — — — — — , —————————— , 


EFQP 

01 


2.0 

1 1. 153SE+7 

0.3 

i _ 

1 • 00(75538 

_ . i 

_ i — — — — 

10,2 

1 

1 


1 3000. 

2.0 

1 

C 

• 

1 

1 

i 

* 

1 

1 


CORD 

13 : 

O 

* 

o 

0.0 

’ 1 

10.0 

1 

t 

11.0 
i _ . 

1 

1 

9 

9 

1 


CORD 

1610.0 

10. 0 

c* 

o 

* 

11.0 

1 

i 

1 

1 

i 


CORD 

14 

10.0 

0. 0 

• 

10.0 

11.0 

_ I 

1 

. _ i _ 

“ i 

* 

i 

« 

< 

i 

iCORD 
1 — — 

15 

10 . 0 

10. 0 

o 

o 

* 

1 1.0 

_ 1 _ 

1 

i 

* 

I 

1 

1 

FIELD DESCRIPTION: 





i - - 

1 


CARD 

FIELD 

FIELD 

CONTENTS 






NO 


LABEL 

NO. 









1 EFQP 
1 NL 

1 E 
1 G 
1 AMU 
1 RHO 
1 ALPHA 

1 AJ 

2 YSTR 
2 NU 

- o CORD 

- 6 IS 

- 6 XENTER 

- 6 YENTER 

- e> ZENTER 

- 6 H 


1 THE ALPHANUMERIC "EPOP H 

1 CODE TO FERFORM LINEAR OR NONLINEAR (ITERATIVE) COMPUTATIONS; 

(INTEGER, 00 * LINEAR, 01 « NONLINEAR) 

3 YOUNG'S MODULUS (REAL > 0.0) 

4 SHEAR MODULUS (REAL > 0.0) 

5 POISSON'S RATIO (REAL, 0.0<AMU<0. 5) 

6 MASS DENSITY (REAL) 

7 MATERIAL PARAMETER IN PLASTICITY LAW (REAL) 

9 MATERIAL PARAMETER IN PLASTICITY LAW (REAL) 

2 NOMINAL YIELD STRESS (REAL) 

3 ALTERNATE POISSON'S RATIO (NOT USED) (REAL) 

1 THE ALPHANUMERIC "CORD" 

1 THE LOCAL GRID NUMBER (INTEGER, 13<IG<16> 

2 THE LOCAL X-COORDINATE OF GRID POINT IG (REAL) 

3 THE LOCAL Y-COGRD I NATE OF GRID POINT IG (REAL) 

4 THE LOCAL Z-COORDINATE OF GRID POINT IG (REAL) 

5 THE THICKNESS OF PLATE AT CORNER CONTAINING LOCAL GRID POINT IG <REAL>0 . ) 


NOTE: 4 CORD CARDS ARE REQUIRED, ONE EACH FOR LOCAL GRID POINTS 13,14,15 AND 16. 

IN THE LOCAL COORDINATE SYSTEM GRID POINT 13 MUST BE AT THE ORIGIN AND POINTS 
13 - 16 MUST ALL LIE IN THE X-Y PLANE, 
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I 


I 


D i 

Z I 

C I X 


I O 1 O I 

I - I . I 

I O I O i 

I -H I -H I 


W I <i I I if) I 

^ | I 1-1 1 r-t J 


•• 

1 

1 

1 

UJ 

1 



1 

! 

H 

1 1 

t 

1 

0 

1 



1 

i 

C 

1 1 

[ 

1 

LL 

1 



1 

f 

z 

1 CL i 

1 Q 

1 

z 

l CL 1 

Q ! 

Q 

! a ; 

1 Q 

X 

1 G I 

1 X 

1 

C 

1 £ ) I 

X ! 

.’r' 

i X I 

1 X 

O 

1 LU i 

i Q 

! 

X 

1 UJ I 

a i 

G 1 

i □ I 

1 O 

Ll 

1 -J 1 

f 0 

i 

LU 

1 0 1 

U ! 

i * i 

w < 

i U i 

1 u 


-J • I 
U J O I 
~ Z I 
U. I 


a _j i 
-J ui i 
uj m i 
H-» <x | 

IL 0 I 


a . i 
x o i 
<r z i 
u i 


- * ’ w ’ W 

Cl O . 
G •* o 

UJ A _| 

-J C A 

. - -J (lJ 
<X llT — f 
cn Lii ^ <x 

0 X LU 

H v Q 

X H-i 

U cn *- 

r d c > 

D Jlth 
ZD ~ 

<r a cn 0 - 


j 

c 

G 

*— » O 
J J J J 

<x c e 

iU Ui Ui CD 

X CL X Z 

^ )-H 

z 

0 CD 3 w 

H H H <1 

H 

fc t t 2 

2 Z Z O 
»— ! <— i »— f LJ 
GOD 
LL CL a. or 
UJ 

a a o z 

w H H 

•X x x a 

0 0 0 U 
Ll Ll II K 

o o o c 


Q UJ 
S iTiTr 

o lu c 
u © z 
= z »~< 
d a 
u z x 
~ a 
x a o 


c x 0 
<10 0 
UJ LU ^ 0 

x x o <r 

H 0 CL Z 


D 0 x 

Z 

C 0 0 

ICC 

ll a u 

DOG 
C 0 0 


LU UJ Lu 
Hhh 
C C C 
Z Z 0 

M W j , 

Q Q 
X X II 
o o o 

0 o 

U U 0 

1 J 0 
>N Li 

Z 

0 0 L: 
•CCU 

u u ►-» 

O D X 

0 0 H 


LU UJ ill LU LU UJ 

X X X X X X 
H H H H- J- h- 


c in >o ^ r*vi ro « 


x x x 

LU LU LU 

CL Q f— f— H- 

a i a ii z z z 

UJ Z X O 0 LU Ui ui 

00 CXU *-* x > nX 


io iii in in in in 


H iH H | 


I f i I 


« m r ; w r - j r * o 

z 


INPUT - LOAD - ELEMENT LOAD DATA 


cj 

ii 


X 

o 


li 

4J 


K 

D 

Cl 

Z 


LU 

CL 

> 

h- 

Z 

o 


c 

D 

G 

LU 

a: 

□ 

u. 


2 

H 

2 


2 

H 

2 


I N i N 


2 

> 

2 

X 


I 

I 

I 

I 2 

i IL 

I Ll 
I Ll! 
I O 
I U 


I > I 
I H I 
! > I 

i 3 r 


> 

> 

> 


A 

> 

2 

x 


2 

Ll 

Ll 

LU 

O 

U 


> 

H 

> 

3 

s 


I > 
I > 
I > 
i x 


! > 

I Ll I 


I LL 


I LU I UJ 
I O I O 
I U I u 


I 


I 3 
I Nl 
I 2 


m i 


~ i 


a i 3 
~ i 


I G 
~ i 


•— I J 
« l 


M 


3 

3 I 

K 

i 

1 


H 

H 1 

D 

i 

1 


— t 
—> 

3 1 

Cl 

i 

1 

> 

SJ 

NI ! 

Z 

i 

1 

w* 

3 

3 1 


i 

1 

0 

> 

> 1 


i 

1 


3 

3 1 

LU 

i 

1 

"W 

X | 

x I 

LL 

i 

1 

iL 

-- - 

- « .. 

> 

— « 

-- 

-- 


1 

h- 

i m 

1 



1 


| M 

1 



1 

LU 

1 c 

1 



1 

-J 

1 *-» 

1 


3 

3 1 


I a 

1 

3 

IL 

Ll 1 

<r ' 

i »— i 

1 


U. 

Ll 1 

K 

1 XI 

1 

0 

lii 

LU 1 


I ►h 

1 

i— i 

G 

G 1 

a: 

1 -M 

1 

w 

U 

1 a i 

o 

1 *“H 

1 

Ll 



1 1 1 

i »■ 

1 

1 

1 

LU 

j 

K 

I 1 1 

1 b- 

1 

I 

1 

-1 

1 

<X 

1 1 1 

1 c 

1 

1 

1 

LL 

1 

r 

1 C 1 1 

i r 

i a 

1 

1 

r 

i a 

a: 

1 C 1 1 

1 X 

1 <r 

i 

1 

<r 

i <X 

o 

i a i i 

1 o 

1 G 

I 

I 

X 

1 O 

Ll 

i -j i i 

! iL 

1 -1 

I 

l 

UJ 

1 -J 


O 


O 

S 

5 


V 

y 

5 


•w‘ 

o 

o 

o 


O 




o 

o 

o 

o 

CJ 


o 

■ 

o 


r*j 
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D 


3 3 
CL O 
Z G 
~ -J 

c 


CD 

H 

Z 

UJ 

K 

Z 

D 

G 



O I 

G • I 
X O f 

H 

.. 1 

~ Z 1 
X I 

1 
\ 


Z 1 

i ~ 


o » 

i 


«-« 1 

Q G 1 


K 1 

G X 1 


X ! 

X X i 


*-« I 

^ <1 1 


X 1 

X G 1 

o 

U 1 

1 

c 

0 1 

1 

o - 

X 1 

/N I 

I 

G *-* 

*-i 1 

i 

I 


i 

O 1 

o . ! 


G » 

X C J 


X ! 

<1 Z 1 


*-< i 

G 1 


X j 

1 




0 


X X X 
X G 

>- h~ 


3 


H <L X 


o 


I G 


G 


Z 1 - G 


G 


o <r 


O 


w LD H 


X 


H X 

<r z x 


1 - 


D H I 


3 


G G f— 


X 


X 


z 

rM 

►“ z 

to 

«— f 


X 3 •-< 



X 

D X 

X 0 

X 

D 

X z z 

o *- 

G 


*-i X 

z 

X 


3 > 

CM X 

D 


O X 

* Z 

X 

*» 

GOG 

— o 


X 

G 

X 

X 

0 X 

O X X 

»» x 

D 

X* 0 

XXX 

X X 


G X 

G <1 

X 

X 

G H 

H Z 

0 O 

K 

<C Z 

- C 3 C 

X z 

a 

f— 

M I 2 h 

K C 

»v 

’W- 

H 

z 


X 

X XX 

^ CO 

O 

O 0 

G 0 I O 

w f— 

z 

z 

X f- 

Z 0 c 

^ 1.1 1 . f 

0 o 

K Z 0 


to 

x 

f—t 
I— J 

CM 


o 

c 

G 

G 


Z ^ 

Li G 

O X 
K L. 

! • i 

rM c 
u 
co - 

LU 

Z !— 

i-i ^ 

-I Cl 
Z 

z « 


-J UJ 


CL 

o 

X 

x 

D 


> 

-J 

Z 


o 


3 

r 


JJD 


> 

-f UJ 
H U 
G X 

Ll’ O 

X Ll 


c 

3 
O 
CO ill 
LU 


a 

<L 

D Q 
€ 


►H C > 


Lii 


J 


Z O D Z 


CO 

rj w 

i— < 

X Ll! 
Ll' G 

r o 

3 G 
z 

C UJ 
X ll 
LL > 

G h- 

c 

uj 3 

X o 
H G 


H IL h Z 
3 > 3 O 
X H X 
- Z Ll 


L 3 w 


!i H !! 

cm to 


o 

0 

Z H 
i-s C 

a o 

*-f G 

a: lu 

G 

0 If 
UJ 

q cm 

Ll' X 

a d 

G 

u 


cn 0 
u z 


to 


K 

<r 

G 

i— < 

O 

z 


CM 

x 

o 


a* 

<r 

H 

a: 

G 

CO 

z 

o 

H 

c 

3 

G 

LU 


UJ 

H 

<E 

a: 


3 

X 

Z 0 

o z 

u a 

i—i 

Z h- 

o <r 

•-» f- 
t- 3 
G X 

z z 

3 O 
IL G 

U Z 

z <r 


3 3 
•-» X 

~ z 

I— < 

0 

Z X 

c o 


0 


<r 

> 

x 

G 

G 

►h 

0 

0 


H X □ 

U UJ X 
Z X _ 
3 Z L! U 
X 3 Z X 
Z ►"* 


z 

X 


UJ 


X 

o 


Q ~ Z 


X X 
G G 
X X 
C G 
X X 

\\ I! 


0 


UJ H 


!! 


CM 


Q 
G 
G ~ 


CM 


H N LO ^ 
Q 
O 
U 


CM 


0 X •-< 
G X 

<r g x 
3 Z X 
G 3 > 
Z X Z H 
X ^ 

X — * X 
X 3 ~ G 
C G 
I! w h C 
* ~ H 
CM ~ 

cox 

^ Z G 
~ C X 


CM 


T"' 


X 

o 

0 

X 

3 

G 

C 

> 


rO 


l ... 


CM 


C 

X 

•—I 

G 

X 

X 

<1 

0 


a . 

C 3 
O G 
G G 
X 
Z ffi 

o 

« X 
H x 
G CD 
c c 

X H 


X 

o 

X 

X 

3 

0 

0 

X 


X 

X 

X 

X 

X 




(c a; 
a. v ) 

H 

1 

i 

« i 

i « 

K 3 

1 

1 i 

i i 

C X 

1 V 

1 ^ 

i 

f 1 
i i 

X X 
H X 

t t 


K 

! I 

l 

o i i 

X O 

1 

f i 

0 

| h 

j 

1 ! 
1 T 

Q « 


1 1 

X O 

<r 

T 

1 

1 CM 

f ^ 

1 I 

i i 

x c 

i 1 
C ! J 

O G 

I 

1 f 

i 

w- | | 

CM 1 1 

X X 
U X 

i 

f ^ 

Vi J I 

I 1 

1 I 

■ <L > 
X t— 

X 

3 X 

1 

t i 

i ' 1 

20 

I i 

i i 
G I 1 

0 G 
X 
X o 

t 1 

I I J 
1 1 

w* j j 

CD ! j 

XX 

H 

1 1 

CM ( 1 

> 

I H | 

-1 .. -G- 

I 1 

0 G 

1 i 

t 1 

i i 

T i 

t i 

i i 

2 . □ 
^ G 
G 

f 1 

! CM 1 

| r- f 

i **< i 

i i 

i i 

« X 

1 i 

n i i 

i i 

X o 
G X 
0 

1 ! 

X o 

« ii 

• — i— j—j. 

Q 

X 



X X 
a t- 



o z 



G X 


CM 


50 


X 

G 

Q 

^ cm r n < 

X 

tn 

LU UJ UJ l 

o X 

131 DOC 

C X 

Q 2 a c 

G X 

UJ UJ UJ L 

X Z 

3 

CD O' O h 

Z Z 

^ ^ C'j 

o 

r r> r 

l-H Q 

N ffi ?> c 

H *H 

^ f *M ^ f* 

O X 

ir* *k i 

<X G 

TM f O *5“ ^ 

X 

► r i 

i— 

^ CM C 

X 

u 

<z 

X 

<-s 

X 

3 cn 

z 

-V 0 

cn x 

X fr- 

X 

O }— 

a x 

H O 

c z 

- X . . 

O 3 

v <r *x 

X Z 

CM . . 

x a 

^ 3 Z Z 

f*M 

<T r-< 

O 

3 X 

r- O 

CD G 

H 

Uj 


X 

N 

X 

rt 

X 

*1 

^ C'M f 1 ; <T 


3 

X 

Z 


> 

fr- 

Z 

Q 


3 

Q 

Ui 

X 

O 

X 


+ 

X! 

I-H 

* 

c 


o 

fr— 

cm 

UJ 

z 


« I 
Z I 
Q I 

i-« I 
fr- I 
X I 
~ I 
X fr 
U I 
Ui I 
X I 
Q 1 
I 

a i 

_j i 


tn 

H 

Z 

UJ 


a 

u 


£~ 

o 

o 

0 

l-H 

H 

l-H 

fr- 

fr- 

fr- 

CI 

ee 

ee 

fr— 

H 

K 

3 

3 

3 

X 

X 

X 

z 

z 

z 

- o 

~ o 

- o 

3 G 

3 U 

3 U 

<r 

ee 

nL 

1 X X 

X X 

X X 

( X 3 

X o 

X o 

t J- 

j 

— X 

«-*» 

1 X z 

X z 

X z 

1 CJ X 

CJ X 

U X 

\ X X 

X X 

X X 

1 O fr- 

O H- 

G H 

1 X 

X 

X 

i ui 

Ui 

Ui 


z 

□ 

1-4 

fr— 

u 

LU 

X 


Z ~ 

o X 

-t H 

H 

a u. 

LU O 
X 


Z *h 
G X 

»-t K 
fr— 

U X 
UJ o 
X 


♦-H ^ 

l-H •* 

l-H r> 

Q > 

CJ > 

Q > 

3 

3 

3 

3 X 

> X 

2 X 

> 

> 

> 

X *-* 

X t-i 

X *-* 

a h 

a fr- 

G fr- 

3 

u 

U 

Z X 

Z X 

z x 

C X 

O X 

a x 

hh cn 

cn 

i-H cn 

H X 

fr- X 

fr- X 

<r x 

<r x 

ee x 

fr— 

— * fr— ++> 

fr— «"> 

3 cn 

O 3 Ui O 

3 U) O 

X X 

X X 

X X 

Z 3 

A Z 3 A 

Z 3 A 

3 uj 

a as 

G 3 


UJ 

tn 

a 

x 


x 

c 


cn 

cn 
z 
o 
►1 

fr- 
<E 
I- 
3 
CL 
Z 

o 
u 

UJ 

o 

X 
G 
X 

X 

a 

LU 


u^tnueetnueecn 
x 


x x lu x x 


u 

X 

z 

X 

UJ 


X X > 


> LU 

•* U 

x X 

o 

LU x 
X 

fr- z 


Z I- z 
X ** *-» 
LU N ^ 
fr— - 
> LU 

cn *cj 

X X 
X o 
H III li 


X 

a 


o 

Ll i-i 
C fr- 

u 

cn uj _ 

i- x *-« 


I- z 
o 

X *+ 
□ fr- 
z CJ 
uj cn uj 
fr- X 


U Z H J 2 H 

M LLj Q M 


o 

u 


3 x 
UJ 

x a 


x* 

z 

D 

X 

x 


a x 


xxx 

O C 3 
X I> X 
fr- 

Z fr- Z 
X •* HH 
X rs» w 
I- •* 

> X 

cn - u 
x X 
X O 
I- X X 

X 

X I- z 
o o 
x ~ 
H O H 
Z CJ 

x cn x 

*-* fr- X 

u z •-< 

w|UO 

it 5 2 

X X 
O x X 
U X □ 


3 • 

i 




X □ 

l-H Z 

i 

t 

C M V) 

m 

<f N 

X 

i 

i 





l 

i 

3 

> 

2 


i 

fr- 

H 

fr- 


fr 

**■ 

•» 

*>■ 

a 3 

i 

3 

> 

2 

3 X 

! 

M 

M 

NJ 

X 3 

i 

3 - 

> ^ 

2 ^ 

~ ee 

1 

X 3 

X > 

X 2 

X 3 

1 

X > 

X > 

X > 


fr 

X - 

X 

X 


fr 

a 3 

□ > 

O 2 


1 

X X 

U X 

U x 


1 

— H T-I 

— H t-H 

— H —H 


1 

+ + 

+ + 

+ + 

9 * 

1 

c c 

c c 

3 C 

X 3 

1 

l-H l-H 

l-H }—f 

l-H l-H 

ee z 

i 

i 

f 1 

1 1 

i j 

a 

1 





I 

M CM 

CM CM 

CM CM 


<x 

X 


cm 


O 

K 


X 

3 



X 


on 

z 


fr- 

»-H 


z 



X 

X 


z 

3 


o 

3 


X 

<E 


X 

fr— 


X 




X 


CJ . 

o 


Z X 

X 


<x > 



o 



cn 3 



fr- <x 

+ 


T 

X 

Ui 

X Q 

z 

fr- 

i-H LU 

l-H 

z 

U 3 

fr- 

X 

l-H l-H 
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2 = PRINT ALL GRID POINT VALUES FOR THIS TIME 

CRITR 3 THE CONVERGENCE CRITERIA FOR ALL DISPLACEMENT VALUES; RELATIVE CONVERGENCE 

DEFAULT TO 0.01 (REAL > O. ) 

DELTM 4 PRINTED OUTPUT TIME INCREMENT 


INPUT - END (REQUIRED LAST INPUT) 
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6.0 OUTPUT OVERVIEW 

Typical output is illustrated in Section 7, Figure 5. An 
echo print of the user input is always printed first. The echo 
print is followed by a group of program controlled output which 
conta ins: 

1. Computed local coordinates of analysis element 

2. Labeled material properties, as input 

3. Tables giving element load values in the local coordi- 
nate system for each input time slice 

4. Table of input initial stress and displacement values. 
The third data output grouping is output each iteration and 
contains the displacement/stress data computed that iteration. 
This output is only at the input time slices. This data is not 
present for the LEQP case since only one iteration is performed. 

The final output consits of displacement/stress data at the 
requested print time increments upon convergence of the problem. 
The printed format is operator controlled for this grouping. 
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7.0 OUTPUT EXAMPLE 

A quadrilateral plate of constant thickness is sued to 
illustrate the input-output format of SFENES. Table 1 provides 
the required input to solve a linear quad plate problem and 
Figure 4 shows the geometric layout of the problem plate. Used 
input data is illustrated in Figure 3 and each input record is 
described in Section 5 of this report. 

Problem output is illustrated in Figure 5. A description of 
this output is contained in Section 6 of this report.. The output 
contrived here is only partial (unbounded D.O.F. only) in order 
to conserve space. 
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TABLE 1 


QUADRILATERAL PLATE MATERIAL PROPERTIES, COORDINATE INPUT, 
LOAD EQUIPMENT AND INITIAL VALUES 

(1) Material Properties 
G * 1 .1 538 X 10 5 * 7 

p - .2588 X 10 -3 
y B .30 

(2) Spatal Coordinates 


Grid No. 

X 

Y 

Z 

Thickness 

13 

0 

0 

0 

1.0 

14 

10.0 

10.0 

0 

1.0 

15 

10.0 

10.0 

0 

1.0 

16 

0.0 

10.0 

0 

1.0 

Time Coordinates 





0.0 

0.0005 

0.0010 


(4) Element Load 

Pressure Load on Top Surface of 

-10,000t or -10,000t in W dir 

(5) Initial Values 

Stress * 0.0 at all grids 

Displacement = 0.0 at all grids 
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500000 IE-03 0.8377444E*02 0.1954736E+03 0. 8377449E+02 0 .2499564E+02 0.2578418E«-02 

9999995E-03 0.1 335305E+03 0.31 15708E+03 0. 1335306E*03 0 . 396429 IE*02 0.6632878£*02 







5000001E-03 0* 19661 OWE *02 0. 197279 1E»C3 0.6657388E+O2 0.7730066E+02 -0.3188586E*02 -0. 9735376E »0l 
9999595G-03 0.3474673E+02 0.3561958E+03 0. H93960E+03 0. 13101318*03 -0.6523181E+02 -0.235581 7E*02 












5000001E-03 0.19 73320E + C3 0.2006250E*02 0.668*766E+02 0. 7686929E+02 -0.91 6853*E*0l -0.3231160€«'02 
O999999E-03 0 .3<»45352E'*'03 0.3594531E+02 0. U99766E*03 0. 13073*IE*03 -0.2265033E+02 -0.6647769E+02 
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5000001E-03 -0.5842536E*00 -0.51003496*00 -0.28892986*01 -0.19326316*00 0.64052966*02 0. 62 76367E *02 
9999999E-03 -0.47902 91E-01 0.16694076*00 0. 5405701E *00 -0.12261966*00 0.98601586*02 0.96912116*02 















5000001E-03 0.0 0.0 0.0 

9999999E-03 0.0 0.0 0.0 
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5000C016-03 0.9298544E+02 0.1378956E+03 0.58063716+02 -0.88449546+02 0. 50430426+02 -0.42756356+02 
99999996-03 0.13566756+03 0.22944496+03 0.98147376+02 -0.15601666+03 0.75294346+02 -0.70683116+02 
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